GET P, N 



= 0 



i 



|0 



1 = 1 + 1 

G = GCD ( P, THOUSAND.PRIMES [i]); 



GETP, M.E 



DETERMINE M-^ MOD P 
BY CALCULATING ' 
GCD (P. M) 



1(0 



RETURN FALSE 



i 



APPLY SIGNED DIGIT 
ALGORITHM TO E 



i 



|6 







RETURN TRUE 





CALCULATE MOD P 



GETU.V 



U3 = U; V3 = V; 
U2 = 0;V2 = 1: 



^'1 



V3 = 0 



T 





RETURN 




US, U2 



Q = LU3/V3j; 
T3 = U3-Q*V3 
T2 = U2-Q*V2 
U3 = V3: V3 = T3 
U2 = V2; V2 = T2 



GETU,V&M 



U3 = U;V3 = V; 



D = 8 4-<M = 32>— ► 
-p 1 Y N /-^ 



D = 5 



[A,B] = TOP2M(U3,V3, M); 




L0G(U3) - L0G(V3) > D 




Q1 = La/bJ; 
T3 = U3-Q1 *V3 



K 




CALCULATE 
XA + YB<2" 




r A 


T3 = XU3 + YV3; 
U3 = V3;V3 = |T3| 



[A. B] = TOP2M(|T3|, V3. M); 
Q2 = Q1- 1a/bJ+1}; 



■<X( [t3|>V3^ 



^3^ 



Q2 = Q1-1 



T3 = U3 - Q2 * V3 



ri 



Q = LU3A/3j 
T3 = U3-Q*V3 



U3 = V3;V3 = T3; 



Y = 1 



Y = GCD (U3,UM0D U3); 
Y = GCD(Y,VMODY); 



RETURN Y 



GETA.B,M . 


1 


r 




M = 32 




Y 



[X, Y] = MOD_EUCLID_GCD(A, B) 



[X, Y] = MOD_LEHMER_GCD(A, B) 



RETURN X.Y 



GET U3, V3, M 



PB1 = LOG (U3) 



BP2 = BP1 - 2M + 1 



RETURN U3[BP1:BP2], V3[BP1:BP2] 



GETU. V 



U3 = U;V3 = V; 
U2 = 0;V2 = 1; 



lj?IO 



Y 




L0G(U3) - LOG (V3) > 32 




[SU,SV] = TOP2M(U3,V3, 16); 
D = LOG (U3)-32; 
J = 0; DONE = FALSE; 
X[0] = 1;X[1] = 0;Y[0] = 0,Y[1] = 1; 



Q = LSU/SVj; 
X[J+2] = X[J]-X[J+1]; 
Y[J+2] = Y[J]-Y[J+1]; 
T = SU-Q*SV;SU = SV;SV = T; 
J=J+1; 



D0NE = SV<-Y[J+1]0R 
SU-SV<X[J+1]-X[J] OR 
SU<2''; 



J EVEN ^ ^^>Y^ 



/pi? 



D0NE = SU<-X[J+1]0R 
SU-SV<Y[J+1]-Y[J]0R 
SU<2°; 



T 



DONE = TRUE?> 



T3 = X[J-1] * U3 + Y[J-1] * V3; T2 = X[J-1] * U2 + YO-1] * V2; 

V3 = X[J] * U3 + Y[J] * V3; V2 = X[J] * U2 + Y[J] * V2; 
U3 = T3;U2 = T2; 



Q = LU3/V3j; 
T3 = U3-Q*V3;T2 = U2-Q*V2 
U3 = V3: V3 = T3; U2 = V2; V2 = T2 

<U3<2^ 



U1 = (U3-U2*V)/U; 
RETURN U1,U2 



GET U, V 



7 f O 



i 



U3 = U;V3 = V; 
U2 = 0;V2 = 1; 



<( jJ3<2^6^ 

i 



7'*/ 



7/6 



U1={U3-U2*V)/U; 
RETURN U1.U2 



Q 


= LU3/V3j; 


T3 = 


U3-Q*V3; 


T2 = 


U2-Q*V2; 


U3 = 


V3;V3 = T3; 


U2 = 


V2;V2 = T2; 



-7(6 



GETU,V&M 



U3 = U;V3 = V; ^<r>j'p 
U2 = 0;V2 = 1; 



[A, B] = TOP2M(U3, V3, M); 




L0G(U3) - L0G(V3) > D 



Y 






Q1 = LA/BJ; 


1 — ► 


T3 = U3-Q1 *V3 







1 




CALCULATE 
XA + YB < 2M 



y 



Q = Q1 



>V3 




[A, B] = TOP2M(|T3|,V3. M); 
Q2 = Q1- (La/BJ + 1) 



Q2 = Q1-1 



T3 = U3-Q2*V3 



T3 = XU3 + YV3; 
T2 = XU2 + YV2: 



T2 = -T2 



Q = Q2 



841 



^#^kT3<0>^ 



U3 = V3;V3 = |T31 
U2 = V2;V2 = T2 



Q = LU3/V3j 
T3 = U3-Q*V3 




T2 


= U2-Q*V2; 




U3 


= V3;V3 = T3; 




U2 


= V2;V2 = T2; 





-^^K^V3 = 0)> 



Y=1 



Y = GCD (U3, U MOD U3); X = LUSA'J; 

COMPUTE X-1 WRT U BY GCD (U, X)- 
U2 = (U2*X-1)M0DU 
Y = GCD(Y,VMODY); 



RETURN Y, U2 



